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1.  Introduction 


In  recent  years,  the  U.S.  Army  Research  Laboratory  (ARL)  researchers  have  achieved  significant 
successes  in  developing  theoretical  models  capable  of  predicting  performance  properties,  such  as 
heats  of  formation  ( 1 ,  2)  and  densities  ( 3 )  of  energetic  materials  (EM).  These  models  allow 
synthetic  and  formulation  chemists  to  safely  and  quickly  screen  candidate  materials  in  order  to 
focus  efforts  only  on  the  most  promising  compounds. 

ARL  researchers  Edward  F.C.  Byrd  and  Betsy  M.  Rice  have  established  correlations  of  quantum 
mechanically  derived  properties  with  heats  of  formation  and  crystalline  densities  for  both  neutral 
and  ionic  solid  energetic  materials.  These  two  properties  are  critical  elements  in  predicting  the 
detonation  and  propulsion  properties  of  a  material.  As  the  crystal  structure  of  a  novel  material 
may  not  be  known,  each  of  the  correlations  requires  quantum  mechanical  data  solely  from  single 
molecules  and  not  any  from  bulk  material.  However,  the  current  codes  used  to  produce  the 
required  information  are  cumbersome  to  use  and  prone  to  failure  unless  strict  formatting  rules  are 
followed. 

In  this  report,  a  brief  summary  of  the  methodologies  employed  will  be  presented,  followed  by  the 
current  state  of  the  codes  used  to  generate  the  final  results.  The  final  goals  of  the  recently 
developed  tools  designed  to  eliminate  user  error  will  be  delineated  and  subsequently  expanded 
upon  in  further  sections.  Finally,  sample  inputs  and  outputs  will  be  presented. 

1.1  Neutral  Heat  of  Formation 

Initially  developed  in  2006  ( /),  the  neutral  heat  of  formation  prediction  tool  uses  Hess’s  law 
(equation  1)  ( 4 )  to  combine  an  atomic  or  group  equivalency  equation  for  the  gas  phase  heat  of 
formation  (equation  2)  coupled  with  an  equation  correlating  properties  of  the  electronic  density 
for  the  heats  of  vaporization  (equation  3)  or  sublimation  (equation  4). 


AH  f  (liquid/sold)  AH  f  (g)  AHvap/sub 

(1) 

AH  f  (g)  =  E  -  Z  ti  f  j 

(2) 

AHvap  =  a(  SAf  +  b(<32,olv) Vl  +  c 

(3) 

AHsub  =  a(S  A)2  +  b(a2totv) 1/2  +  c 

(4) 

In  equation  2,  E  is  the  quantum  mechanically  determined  electronic  energy,  nt  is  either  the 
number  of  atom  or  group  types  j  contained  in  the  neutral  molecule,  and  sj  is  either  the  atom  or 
group  equivalent  energy  of  atom/group  j.  These  equivalent  energies  were  determined  through 
the  fitting  of  experiment  gas  phase  heats  of  formation  [AH°f  (g)]  to  quantum  mechanically 
computed  energies  E. 
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Equations  3  and  4  use  statistics  derived  from  the  electrostatic  surface  potential  (ESP),  previously 
demonstrated  by  Politzer  et  al.  (5-7)  to  correlate  to  heats  of  vaporization  and  sublimation. 
Following  the  conventions  lay  down  by  Politzer  et  al.  (8)  and  Rice  et  al.  (9),  the  definitions  of 
these  electrostatic  properties  are  global  properties  across  the  0.001  electron/bohr  isosurface  of 
electronic  density.  The  electrostatic  potential  V  is  defined  as: 


V(r)=Y 

^  J  ZjjlR;-rl  J  Ir'-rl  ’ 


(5) 


where  Z,  and  R,  represent  the  charge  and  position  of  atom  i,  and  p( r)  is  the  electronic  density. 
The  statistical  tools  developed  by  Politzer  et  al.  (8,  10)  used  in  equations  3  and  4  are  the  surface 
area  (SA),  the  variance  of  electrostatic  surface  potential  (o'),  and  the  balance  parameter  (v). 
These  are  defined  as: 


erv 


7  _ 2  i 

°tot  —  °+  +  cr_ 


v  — 


KO2 


(6) 

(7) 

(8) 
(9) 


where  E(r,j  is  the  value  of  the  ESP  at  r,  and  X  is  the  average  electrostatic  surface  potential  (x  is 
either  the  positive  (+)  or  negative  (-)  portion  of  the  ESP).  The  surface  area  is  the  surface  area  of 
the  0.001  electron/bohr3  isosurface  of  the  electron  density  of  the  molecule,  a2tot  is  a  measure  of 
the  variability  of  electronic  potential  on  the  surface,  and  v  is  the  degree  of  balance  between  the 
positive  and  negative  charges  on  the  isosurface.  Equations  3  and  4  each  come  with  a  different 
set  of  fitted  a,  b,  and  c  parameters. 


As  illustrated  in  reference  1,  this  methodology  enjoys  considerable  accuracy,  with  the  root  mean 
square  error  (rms)  for  the  solid  phase  heat  of  formation  under  6  kcal/mol  (figure  1).  Table  1 
reveals  the  excellent  agreement  with  experimental  gas-,  liquid-,  and  solid-phase  heats  of 
formation  this  methodology  exhibits,  with  the  group  equivalent  energies  yielding  slightly  better 
solid-phase  heats  of  formation  over  the  atom  equivalent. 
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Experimental  Values  (kcalmol) 


Figure  1.  Accuracy  of  the  neutral  solid  phase  heat  of  formation  tool  (reference  1). 

Table  1.  Root  mean  square  errors  and  maximum  deviations  for  the  gas, 

liquid  and  solid  phase  heat  of  formation  for  neutral  materials  using 
the  B3LYP/6-31 1++G(2df,2p)  energies  (reference  1). 


Atom  Equivalent 

Group  Equivalent 

rms 

error 

max  error 

rms  error 

max  error 

AH°f(£) 

2.8 

5.4 

3.2 

6.5 

AH°fm 

3.1 

7.2 

3.2 

7.4 

6.0 

16.7 

5.6 

12.2 

The  current  methodology  required  to  compute  the  neutral  solid-phase  heat  of  formation  for  a 
notional  material  is  as  follows: 

1.  Build  the  notional  molecule. 

2.  Compute  the  optimal  geometry  using  the  B3LYP  spin-restricted  Kohn-Sham  density 
functional  theory  (KS-DFT)  (11)  with  the  6-31G*  Pople  gaussian  basis  set  (12),  currently 
done  using  the  Gaussian09  program  package  (13). 

3.  Check  whether  the  final  geometry  is  a  stable  minimum  or  if  it  is  a  transition  state 
(determine  from  an  additional  Gaussian09  vibrational  frequency  calculation). 

4.  Compute  the  electronic  energy  of  the  optimized  geometry  (step  2)  using  the  B3LYP  density 
functional  theory  with  the  6-311++G(2df,2p)  Pople  gaussian  basis  set  (12). 
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5.  Build  the  cube  and  pot  files  from  the  wavefunction  data  using  the  Gaussian  formchk  tool. 

6.  Run  the  ARL  ED  AT  tool,  developed  for  ARL  researchers  by  Jennifer  Hare,  to  determine 
the  molecular  surface  area  along  with  the  various  required  Politzer  parameters. 

7.  Determine  the  number  of  atom/group  types  for  the  notional  molecule. 

8.  Plug  the  aforementioned  computed  parameters  into  the  neutral  heat  of  formation  formula. 

This  is  a  fairly  lengthy  process  and  is  prone  to  transcription  errors.  However,  all  this  is  required 
to  generate  just  one  of  the  potentially  desired  properties,  and  this  procedure  is  not  the  most 
complicated  presented  herein. 


1.2  Crystalline  Density 

Developed  in  2007  (5),  the  crystalline  density  predictor  methodology  correlates  the  quantum 
mechanically  derived  electronic  density  volume  of  the  molecular  unit  to  the  crystalline  density. 
Similarly  to  the  neutral  heat  of  formation  correlation,  this  correlation  also  has  exhibited 
significant  success  in  predicting  the  crystal  density.  As  discussed  in  reference  3,  it  was 
determined  that  the  crystal  density  correlated  very  well  with  the  quantum  mechanically  derived 
molecular  volume.  Rice  et  al.  found  that  the  0.001  electron/bohr  isosurface  of  the  electron 
density  of  the  isolated  molecule  (the  same  isosurface  as  used  for  the  neutral  heat  of  formation 
tool)  encapsulated  the  molecular  volume  and  yielded  an  excellent  correlation  to  density  (figure 
2). 


Figure  2.  Accuracy  of  the  crystalline  density  prediction  tool.  The  image  on  the  left  is 

for  neutral  species,  while  the  image  on  the  right  is  for  ionic  compounds.  Associated 
root  mean  square  errors  are  presented  (reference  3). 

The  derived  correlation  contains  no  fitted  parameters,  and  simply  requires  the  quantum 
mechanically  derived  molecular  volume  and  the  molecular  mass  (equation  10). 

Pcrystai  =  (total  molecular  mass)  /  (total  molecular  volume)  (10) 
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While  it  is  simple  to  determine  the  molecular  volume  for  neutral  species  (as  there  is  only  one 
molecule),  calculating  the  formula  unit  volumes  of  ionic  systems  can  be  difficult  if  no  a  priori 
knowledge  of  the  crystal  structure  exists.  As  a  multi-component  system  (e.g.,  ionic  systems)  can 
have  a  near  infinite  combination  of  potential  relative  separations  and  orientations  of  the 
components  in  the  formula  unit,  it  can  be  problematic  to  determine  the  most  probable 
arrangement  that  most  closely  corresponds  to  the  crystalline  phase  without  this  advanced 
knowledge.  Any  attempt  to  sample  the  phase  space  of  possible  configurations  has  no  guarantee 
of  properly  capturing  the  correct  experimental  crystal  structure.  In  order  to  avoid  testing 
multiple  different  configurations,  one  can  assume  additivity  of  the  individual  components  to 
define  formula  unit  volumes  (14-16). 

Assuming  an  additive  approach  implies  that  the  volume  of  the  formula  unit  MpX9  of  an  ionic 
crystal  is  simply  the  sum  of  the  volumes  of  the  ions  contained  in  the  formula  unit: 

Volume  =  pVM+  +  gVX-,  (11) 

with  M  denoting  the  cation  and  X  the  anion.  For  formula  units  containing  more  than  two 
components  (such  as  when  solvent  molecules  are  present  in  the  crystal  lattice),  the  additional 
molecules  are  also  included  in  the  total  volume: 

V  olume  =  pW  M+  +  <r/V X-  +  rVsoivent ,  (12) 

where  r  is  the  number  of  solvent  (or  additional)  molecules  in  the  formula  unit. 

Rice  et  al.  determined  that  for  ionic  systems,  errors  were  introduced  with  increasing  hydrogen 
content,  such  that  there  was  a  linear  correlation  with  the  number  of  hydrogen  atoms  in  an  ion  and 
the  error  in  the  predicted  volume.  It  was  concluded  this  error  was  due  to  overly  long  X-H  (X  = 
N,  C)  bonds,  as  predicted  by  the  quantum  mechanical  method  employed.  Therefore,  a  correction 
to  the  molecular  volume  was  established  in  order  to  reduce  the  error  in  the  initial  quantum 
mechanical  volume: 


V corrected  —  Vunc0rrected  —  [0.6763  +  0.9418  (#  Hion)]  (13) 

The  current  methodology  for  calculating  the  crystalline  density  for  a  notional  material  is  as 
follows: 

1.  Build  the  notional  molecule. 

2.  Compute  the  optimal  geometry  using  the  B3LYP  KS-DFT  with  the  6-3 1G**  Pople 
gaussian  basis  set  (12). 

3.  Check  whether  the  final  geometry  is  a  stable  minimum  or  if  it  is  a  transition  state 
(determine  from  an  additional  vibrational  frequency  calculation). 

4.  Build  the  cube  and  pot  files  from  the  wavefunction  data  using  the  Gaussian  formchk  tool. 
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5.  Run  the  ARL  ED  AT  tool  to  determine  the  molecular  volume  and  correct  for  number  of 
hydrogens. 

6.  Determine  the  molecular  mass. 

7.  Plug  the  computed  parameters  into  the  density  formula. 

1.3  Ionic  Heat  of  Formation 


The  ionic  heat  of  formation  tool  was  developed  in  2009  (2)  and  followed  a  similar  methodology 
as  that  for  the  neutral  heat  of  formation  correlation.  However,  for  ionic  materials,  one  cannot  use 
the  heat  of  sublimation  coupled  with  the  gas  phase  heat  of  formation  (as  in  equation  1),  as  the 
intermolecular  forces  binding  the  ions  together  are  entirely  different  than  those  forces  found  in 
neutral  materials.  The  methodology  employed,  a  modified  form  of  the  Born-Fajans-Haber  cycle 
(4),  is  formulaically  similar  to  equation  1: 


M(g)  +X(g) 


IE,EA 


■»M(g)  +Xfe) 


AH, 


Lattice 


->M+X 


solid 


(14) 


where  instead  of  performing  the  full  cycle  (equation  14),  the  gas  phase  heats  of  formation  are 
determined  directly  for  the  ions,  without  ever  computing  any  information  on  the  neutral 
precursor  species  (equation  15). 


AH 


Lattice 


-»M+X 


solid 


(15) 


While  the  gas  phase  heats  of  formation  are  computed  quantum  mechanically,  the  lattice  enthalpy 
( AHUtticc  )  is  calculated  from  the  lattice  potential  energy  (Upot),  a  measurement  of  the  binding 

energy  of  the  ionic  crystal  (77): 


AHL^ce  =  Upot(M;Xq )  +  [p(  V/2  -  2)  +  q(nx./2  -  2)]RT 


(16) 


where  nM+  and  nx  depend  on  whether  the  ions  are  monatomic,  linear  polyatomic,  or  non-linear 

polyatomic  species,  and  T  is  298  K.  As  discussed  in  reference  2,  the  calculation  of  the  lattice 
potential  energy  will  use  either  the  Jenkins  model  (14,  18,  19)  or  the  Gutowski  model  (20), 
which  is  applicable  only  for  1:1  salts.  Both  the  Jenkins  and  Gutowski  models  correlate  the 
molecular  volume  to  the  lattice  enthalpy: 


Upo,=  21 


W+P]’ 


(17) 


where  I  is  the  ionic  strength,  V  is  the  volume,  and  a  and  P  are  empirically  derived  parameters 
(these  differ  between  the  Jenkins  and  Gutowski  methods).  Please  note  that  the  hydrogen- 
corrected  molecular  volumes  are  used  when  computing  the  molecular  volume  (equation  13). 
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As  illustrated  in  table  2,  the  methodologies  do  reasonably  well  for  predicting  the  solid-phase 
heats  of  formation  for  ionic  materials.  The  errors  are  larger  than  those  observed  for  the  neutral 
heat  of  formation  tool;  however,  one  must  note  that  the  binding  energies  for  ionic  materials  span 
a  much  larger  range  (-100  kcal/mol)  than  those  usually  seen  for  neutral  materials  (~10  kcal/mol) 
While  marginally  better  results  can  be  obtained  for  non-l:l  salts,  assuming  one  had  knowledge 
of  the  crystal  structure,  the  Jenkins  methodology  allows  for  a  completely  predictive  tool  for  the 
determination  of  solid-phase  heats  of  formation  of  molecular  energetic  salts  without  requiring 
this  information.  Surprisingly,  the  Gutowski  method  for  1:1  salts  yielded  the  best  overall  results 
(2). 


Table  2.  Solid  phase  heats  of  formation  for  ionic  materials  predicted  using 
various  methods  to  calculate  component  AH°f(g)  and  AHLattice 
(reference  2). 


Theoretical  Method 

AH°f(g)  +  AH°Lattice 

unsigned  mean 
error  (kcal/mol) 

rms  error 

(kcal/mol) 

1:1,  2:1,  2:2  salts 

G3MP2(B3LYP)  +  Jenkins 

29.7 

36.6 

1:1  salts 

G3MP2(B3LYP)  +  Gutowski 

19.1 

24.0 

G3MP2(B3LYP)  +  Jenkins 

24.4 

28.0 

The  current  methodology  for  calculating  the  ionic  heat  of  formation  for  a  notional  material  is  as 
follows: 

1.  Build  each  of  the  components  (subunits)  in  the  notional  molecular  system. 

2.  Compute  the  G3MP2(B3LYP)  (21)  gas-phase  energies  for  each  subunit. 

3.  Check  whether  the  final  geometries  are  stable  minima  or  if  they  are  a  transition  state 
(determined  from  the  vibrational  frequency  calculation  included  in  the  G3MP2(B3LYP) 
calculation). 

4.  Compute  the  optimal  geometry  using  the  B3LYP  KS-DFT  with  the  6-3 1G**  Pople 
gaussian  basis  set  for  each  subunit. 

5.  Check  whether  the  geometries  computed  with  B3LYP/6-31G**  are  stable  minima  or  if 
they  are  transition  states  (determine  from  an  additional  vibrational  frequency  calculation). 

6.  Build  the  cube  and  pot  files  from  all  the  wavefunction  data  using  the  Gaussian  formchk 
tool. 

7.  Run  the  ARL  EDAT  tool  to  determine  the  molecular  volume  of  each  of  the  B3LYP/6- 
31G**  structures  and  correct  for  hydrogen. 

8.  Compute  the  lattice  potential  from  the  total  system  volume  using  either  the  Jenkins  or 
Gutowski  methodology  (dependent  on  whether  the  system  is  a  1:1  salt). 
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9.  Compute  the  lattice  energy  from  the  lattice  potential. 

10.  Compute  the  total  gas-phase  heat  of  formation  from  the  G3MP2(B3LYP)  values. 

11.  Calculate  final  total  solid-phase  heat  of  formation  from  gas-phase  energies  and  lattice 
energy. 

1.4  Design  Goal 

As  illustrated  in  the  previous  sections,  the  current  state  of  the  tools  requires  an  elaborate  series  of 
calculations  to  be  performed  in  a  specific  order.  Unmentioned  are  the  idiosyncrasies  inherent  in 
the  in-house  codes  used  to  compute  the  required  information  from  the  Gaussian  output  files. 
There  are  several  issues  with  the  current  toolset: 

1.  Tools  are  not  coupled;  i.e.,  must  perform  distinct  calculations  for  each  property  separately. 

2.  Not  “push  button”  ready;  i.e.,  significant  effort  required  prior  to  even  generating  necessary 
initial  Gaussian  data. 

3.  Require  extremely  format  specific  inputs.  For  example,  an  ARL  researcher  discovered  that 
being  off  by  a  single  space  or  carriage  return  caused  immediate  failure. 

4.  In  general,  the  overall  process  is  tedious  and  rife  for  failure,  especially  for  the  non-expert 
user  (as  can  be  expected  from  a  series  of  codes  built  up  piecemeal  in  order  to  solve  specific 
issues). 

The  final  design  goal  of  any  solution  designed  to  address  each  of  these  issues  would  need  to 
have: 

1.  A  simple  input  structure,  such  as: 

(a)  xy  z  coordinates  and  atom  type  (C,  H,  N,  etc.); 

(b)  a  single  line  of  information  per  molecule  in  the  overall  system  (e.g.,  charge  of  ionic 
partners,  number  of  saturated  carbons,  etc.);  and 

(c)  as  freeform  an  input  as  possible,  limiting  possible  sources  of  failure  due  to  formatting 
issues. 

2.  A  means  of  monitoring  all  required  calculations,  limiting  user  input  to  an  initial  input 
structure,  and  receiving  the  final  computed  results. 

3.  A  coupled  procedure,  where  the  researcher  need  only  submit  one  job. 

Each  of  these  needs  has  successfully  been  addressed  and  shall  be  expounded  upon  in  the 
following  section. 
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2.  Method  of  Solution 


2.1  Code 

Three  bash  scripts  have  been  written  to  address  the  previously  listed  issues: 

1.  One  main  script  that  drives  the  other  scripts,  as  well  as  submits  the  required  Gaussian 
calculations  and  monitors  the  queues  (master_densheat.s) 

2.  A  script  to  build  the  required  Gaussian  input  files  (gbuild.s) 

a.  For  density  calculations,  the  script  will  build  a  file  called  filename jdensity. com  if  the 
molecule  is  a  neutral  or  filename_  [ion]_[number]_density.com  for  ionic  species 
(where  [ion ]=cation,  anion ,  or  neutral  and  [number]  is  the  running  count  of  that  ion 
type  (one  through  n)). 

b.  For  heat  of  formation  calculations,  the  script  will  build  a  file  called  filename  com 
if  the  molecule  is  a  neutral  or  filename_  [ion]_[number]_heat.com  for  ionic  species 
(where  [ion ]=cation,  anion,  or  neutral  and  [number]  is  the  running  count  of  that  ion 
type  (one  through  n)). 

3.  A  script  to  construct  the  data  needed  to  determine  the  heats  of  formation  and/or  crystalline 
densities  and  subsequently  calculating  said  properties  (arl_densheat.s) 

Given  a  properly  formatted  input  file,  the  scripts  will  analyze  the  input  file  and: 

•  Load  the  relevant  modules  (master_densheat.s,  arl_densheat.s); 

•  Determine  appropriate  format  needed  and  build  the  required  Gaussian  input  files  (gbuild.s); 

•  Submit  the  built  input  files  and  monitor  the  PBS  queue  to  determine  when  calculations  are 
complete  (master_densheat.s); 

•  Copy  the  required  Gaussian  output  files  from  the  scratch  directory  to  the  working  directory 
(master_densheat.s); 

•  Generate  the  required  Gaussian  utility  formchk  .cube  and  .pot  files  (arl_densheat.s); 

•  Run  the  ED  AT  tool  in  order  to  generate  required  Politzer  statistics  (arl_densheat.s); 

•  Gather  required  data  from  all  generated  files  (e.g.,  surface  areas,  energies,  volumes,  etc.) 
(arl_densheat.s);  and 

•  Compute  requested  heats  of  formation  or  densities  (arl_densheat.s). 

The  scripts  also  perform  a  limited  series  of  error  checking: 
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•  Checks  that  for  optimized  structures,  all  frequencies  are  real  (i.e.,  the  structure  is  not  a 
transition  state  structure,  but  rather  a  minima);  if  not,  the  script  will  stop  and  will  print  an 
error  message  in  the  arl_heat_density.err  error  file.  (master_densheat.s) 

•  Checks  that  the  required  Gaussian  outputs  do  not  already  exist  (the  .chk  and  .out  files),  and 
if  they  do,  will  not  resubmit  these  Gaussian  calculations.  Note  that  if  modifications  were 
made  to  the  .dat  input  file,  it  is  necessary  to  delete  previous  .chk  and  .out  files  remaining  in 
the  working  directory.  (master_densheat.s) 

•  For  ionic  systems,  scripts  check  the  total  charge  is  neutral  (original  tools  were  not 
developed  for  non-net-neutral  ionic  systems);  if  not,  the  script  will  halt. 

(mas  ter_densheat .  s ) 

•  Checks  that  the  Gaussian  calculations  ran  successfully  to  completion.  An  error  message 
will  be  printed  in  the  arl_heat_density.err  error  file  if  Gaussian  failed  to  complete. 
(arl_densheat.s) 

•  For  neutral  calculations,  the  script  will  check  that  the  number  of  atom  types  listed  in  the 
first  line  (e.g.,  number  of  carbons  equals  C  +  C’  atom  types)  correspond  to  the  number  of 
atoms  in  the  Cartesian  coordinates.  Currently,  there  is  no  error  check  to  determine  if  the 
number  of  group  types  is  correct  (e.g.,  C-N02,0-N02,  etc.)  If  there  is  a  disagreement,  the 
script  will  halt. 

For  large  molecular  systems,  an  additional  master  script  is  available  ( large_master_densheat.s ), 
which  acts  exactly  the  same  as  the  standard  master_densheat.s  script  with  the  notable  exception 
of  requesting  64  cores  and  96  h  of  queue  for  use  by  the  Gaussian09  code.  Even  larger  amounts 
of  CPU  time  and  a  larger  number  of  cores  can  be  requested  through  personal  modification  of  a 
local  copy  of  the  script  or  by  request. 

2.2  Gaining  Access  to  Required  Codes 

The  ARL  scripts  are  located  on  the  ARL  DSRC  machine  “Harold,”  with  the  host  name  used  to 
log  on  is  “harold.arl.hpc.mil.”  The  directory  the  scripts  are  located  in  is  “/usr/people/ebyrd/bin/”, 
which  can  be  aliased  if  so  desired. 

In  order  to  use  these  scripts,  prospective  users  must  have  access  to  the  ARL  DSRC,  as  well  as 
access  to  the  Gaussian  software  package.  Access  to  the  ARL  DSRC  can  be  obtained  through  the 
following  route  (current  as  of  April  2012,  subject  to  change): 

1.  Go  to  the  Portal  to  the  Information  Environment  at 
https  ://ieapp  .erdc  .hpc  .mil/info/kerberos  V  alidate.  j  sp 

2.  At  the  bottom  of  the  grey  box  in  the  middle,  click  on  “Request  login  account”. 

3.  Click  on  “Apply  for  an  IE  Account”  and  fill  out  the  required  information.  The  required 
Kerberos  Realm  is  “HPCMP.HPC.MIL”  for  US  government  employees  or  for  those  who 
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intend  to  get  a  NAC.  For  foreign  nationals  or  others  who  do  not  want  a  NAC,  access  can 
only  be  obtained  for  the  Open  Systems,  and  in  that  case,  the  Kerberos  Realm  would  be 
ARSC.EDU.  Plans  are  in  place  to  transition  these  scripts  to  the  Open  Resource  System. 

4.  Select  the  organization  ID  “ARLAP”. 

5.  The  rest  of  the  questionnaire  is  self-explanatory  and  should  be  completed. 

6.  Next,  go  back  to  https://ieapp.erdc.hpc.mil/info/kerberosValidate.isp 

7.  At  the  top  right-hand  corner,  click  on  “IAA  training.”  Complete  the  training,  and  print  out 
and  sign  the  certificate,  and  either  FAX  or  scan/email  to  the  ARL  S/AAA  (Service/Agency 
Approval  Authorities)  (Dr  Betsy  Rice  is  alternative  S/AAA  and  will  accept  certificates). 

To  obtain  access  to  Gaussian,  contact  the  CCAC  (Consolidated  Customer  Assistance  Center) 
(http://www.ccac.hpc.mil/)  and  request  access  on  the  appropriate  DSRC  (currently  ARL)  by 
emailing  the  help  desk  at  Help@ccac.hpc.mil. 

Prior  to  running  the  scripts  for  the  first  time,  it  is  advised  to  test  that  access  has,  indeed,  been 
granted  for  all  required  codes.  This  is  only  required  as  a  test  to  determine  if  proper  access  has 
been  granted.  It  is  not  required  to  do  this  each  time  prior  to  using  the  script. 

To  test  if  proper  access  has  been  granted,  at  the  prompt  (denoted  “»”  in  the  text)  type  the 
following  commands: 

»module  load  g09_c01 
»which  g09 

where  the  g09_c01  above  is  g-zero-nine-underscore-c-zero-one.  The  which  command  should 
return  something  other  than  “g09:  Command  not  found.”  Subsequently  type  the  following 
commands: 

»module  use  /usr/cta/CSE_MSRM_ATOMISTI C/modules 
»module  load  cseinit-msrm-atomistic  cse-msrm/atomistic/edat/latest 
»which  edat 

where  the  which  command  should  return  something  other  than  “edat:  Command  not  found.”  If 
either  of  the  two  commands  fails,  please  contact  the  CCAC  help  desk. 

2.3  Usage 

The  ARL  scripts  require  only  one  input  file  per  molecular  system  (referred  to  a.%  filename AdX  in 
the  text).  The  file  name  must  end  in  “.dat”  and  contain  no  other  periods  or  special  characters 
aside  from  and  which  are  acceptable.  Multiple  molecules  can  be  submitted  at  one  time. 

»  /usr/people/ebyrd/bin/master_densheat.s  filenamel.dat  filename2.dat  etc. 
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For  simplicity,  in  this  and  all  subsequent  examples,  the  call  for  the  script  is  assumed  either  to  be 
aliased  to  or  explicitly  entering  the  current  path  (i.e.,  /usr/people/ebyrd/bin/master_densheat.s), 
and  shall  be  shortened  in  the  text,  as  illustrated  for  the  previous  example: 

»master_densheat.s  filename l.Astt  Jilename2.dat  etc. 

If  the  script  has  not  been  aliased  to  the  full  path,  or  the  full  path  is  not  used,  the  script  will  not 
function.  Please  contact  CACC  with  aid  in  setting  up  an  alias  or  type  out  the  full  path  for  each 
use  of  the  script  (/usr/people/ebyrd/bin/master_densheat.s). 

The  user  can  define  either  a  density  calculation  or  heat  of  formation  calculation  (default  is  to  do 
both)  using  the  “-d”  and  “-h”  flags: 

»master_densheat.s  -h  filename. dat 
»master_densheat.s  -d  filename. dat 

The  flags  that  may  be  called  are  as  follows: 

-b  Only  build  the  required  Gaussian  input  files. 

-b 

-build 

-build 

-h 
-h 
-heat 
—heat 

-d  Only  compute  the  crystalline  density  and  writes  data  to 

— d  density.out. 

-density 

—density 

-help  Print  help  page  for  script  to  screen. 

The  script  may  be  run  in  the  background  (using  the  &  command),  allowing  continued  use  of  the 
xterm  window  (as  well  as  allowing  the  script  to  run  if  the  xterm  window  is  closed),  with  standard 
out  data  printed  to  a  file  through  the  following  command: 

»master_densheat.s  filename. dat  >  tempname  & 

If  the  background  command  (&)  is  not  used  upon  logging  out  of  the  xterm  window,  the  script 
will  terminate  and  not  produce  the  final  results.  It  is  strongly  recommended  that  the  background 
command  be  used  in  all  submissions  in  order  to  prevent  re-running  of  the  script  due  to  the  xterm 
window  prematurely  closing. 


Only  compute  the  heat  of  formation  and  write  data  to 
heat_of_formation.out. 
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2.4  Format  of  the  Input  File 

The  input  file  (filename. dat)  will  have  one  of  two  formats — one  for  neutral  molecules,  and  one 
for  ionic  molecules.  The  script  will  determine  the  appropriate  files  to  construct  and  submit  based 
on  the  data  in  the  input  file.  The  input  files  are  almost  completely  format-free,  with  the  only 
caveat  that  the  first  line  must  not  be  blank. 

For  neutral  molecules,  the  first  line  must  delineate  the  atom/group  breakdown  of  the  molecule,  as 
per  reference  1.  This  required  data  must  be  input  with  the  following  order: 

CHNOC'N’O'  C-N02  N-NO2  O-NO2  C-N3  X-NO  , 

with  the  first  four  (CHNO)  denoting  the  number  of  unconjugated  carbon,  hydrogen,  nitrogen, 
and  oxygen  atoms  in  the  molecule.  The  following  three  (C’N’O’)  are  the  number  of  conjugated 
carbon,  nitrogen,  and  oxygen  atoms  in  the  molecule.  The  remainder  are  the  group-based 
additivity  values,  with  the  first  three  the  number  of  nitro  groups  (NO2)  bonded  to  carbon, 
nitrogen,  and  oxygen,  followed  by  the  number  of  azido  groups  (N3)  bonded  to  carbon  and  the 
total  number  of  nitroso  moieties  (NO).  Using  nitromethane  (CH3NO2)  as  an  example,  this  would 
yield  as  a  first  line  to  the  input  file: 

130001210000 

Subsequent  lines  in  the  input  file  must  be  atom  type  (either  atomic  number  or  symbol  (e.g.,  C  or 

o 

6)),  followed  by  the  Cartesian  coordinates  in  Angstroms.  The  number  of  spaces  between  the 
atom  types  and  coordinates  is  not  fixed  (free  format).  Additionally,  blank  lines  between  rows 
will  be  ignored. 

Similarly  for  ionic  materials,  the  first  line  must  denote  first  the  charge  of  the  moiety,  followed  by 
the  number  of  times  the  ionic  moiety  appears  in  the  total  molecular  complex.  Subsequent  lines 
must  be  atom  type  (either  atomic  number  or  symbol  [e.g.,  N  or  7]),  followed  by  the  Cartesian 

o 

coordinates  in  Angstroms  (once  again  format  free  and  ignoring  blank  lines).  This  format  is 
repeated  for  each  moiety  in  the  molecular  complex.  The  order  of  the  moieties  within  the  input 
file  is  irrelevant;  however,  the  total  charge  must  be  zero  for  the  molecular  system. 

The  Appendix  gives  a  brief  tutorial  on  the  gaussview  program,  which  can  be  used  to  construct 
Cartesian  coordinates  if  needed. 

2.5  Examples 

Sample  test  cases  with  associated  inputs,  Gaussian  outputs,  and  final  density  and  heat  of 
formation  outputs  are  included  in  the  /usr/people/ebyrd/bin/script_examples/  directory. 
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Neutral  example:  Nitromethane  (CH^NCh): 


130001210000 
C  -1.91299700  -0.32347600 
7  -0.70935800  0.57063300 

O  0.08637700  0.44472600 
O  -0.62778900  1.38491900 
1  -2.63104348  0.12211454 

1  -2.32224168  -0.35240633 

H  -1.60067525  -1.30351263 


7.60644100 

7.59458100 

8.51978100 

6.68070600 

8.29868770 

6.59826669 

7.96292804 


Ionic  example: 

-21 


triaminoguanidinium  bis(triazolyl)amine  [(CH<)Nf,+h  (C2  H  N9)  (H20)2] 


N 

4.688  3.109 

2.826 

7 

4.211  3.610 

1.647 

N 

4.127  4.909 

1.734 

N 

4.533  5.323 

2.956 

N 

5.323  4.200 

4.878 

N 

5.368  1.836 

5.274 

N 

5.855  1.085 

6.297 

N 

6.373  1.863 

7.188 

7 

6.246  3.161 

6.792 

C 

4.855  4.191 

3.586 

'4  « 

6 

5.641  3.069 

5.624 

H 

5.570  4.924 

5.253 

>•  >*** 

02 

*  f 

O 

7.781  -0.916 

3.731 

n  j 

H 

7.872  -1.798 

3.449 

H 

12 

8.121  -0.717 

4.463 

N 

N 

3.176  -0.379 
4.003  0.041 

2.305 

3.372 

••  f 

N 

5.066  -1.276 

1.348 

N 

5.645  -1.951 

0.251 

»w 

< 

1 

m 

N 

3.019  -1.562 

0.333 

N 

1.613  -1.420 

0.427 

)  j*  . 

C 

3.742  -1.085 

1.334 

J9 

H 

1 

3.416  -2.212 
1.273  -2.195 

-0.184 

0.685 

1 

2.196  -0.318 

2.304 

H 

6.301  -2.575 

0.579 

1 

4.119  0.937 

3.383 

H 

5.555  -0.738 

1.869 

H 

3.557  -0.287 

4.134 

H 

1.231  -1.224 

-0.382 

H 

6.086  -1.279 

-0.158 

14 


2.6  Restrictions  and  Errors 


Please  note  that  there  are  some  restrictions  on  the  atom  types  capable  of  being  handled  by  the 
various  tools.  For  neutral  materials,  the  heat  of  formation  tool  has  only  been  parameterized  for 
CHNO  containing  molecules.  For  ionic  materials,  while  the  capability  exists  for  ionic  heat  of 
formation  calculations  for  molecules  containing  hydrogen  (H)  through  to  chlorine  (Cl),  the  tool 
has  only  been  validated  for  CHNO  materials,  and  is,  therefore,  restricted  to  these  atoms.  It  is 
possible  to  expand  the  range  of  atoms  for  the  ionic  heat  of  formation  calculation  by  editing  the 
arl_densheat.s  script,  and  commenting  out  the  current  declare  -a  g03e0  line  near  the  top  of  the 
script  and  commenting  in  the  subsequent  declare  -a  g03e0  call.  Meanwhile,  though  the  density 
predictor  tool  has  no  inherent  limitations  on  atom  type,  this  tool  has  only  been  validated  for 
CHNO-containing  compounds. 


Currently,  the  script  assumes  the  Gaussian  job  will  complete  within  4  h  using  32  processors  and 
a  ptile  of  8.  Additionally,  the  time  between  checking  the  qstat  run  queue  is  set  to  5  min.  If  these 
values  are  insufficient  due  to  the  size  of  the  calculation  required  (most  likely  due  to  large  ionic 
calculations),  then  there  exists  the  large_master_densheat.s  script,  which  requests  96  h  and  64 
CPUs.  If  this  script  is  insufficient,  then  modification  of  the  script  is  required.  The  relevant 
parameters  can  be  set  at  the  top  of  a  local  copy  of  the  master _densheat.s  script.  The  following 
parameters  can  be  modified  at  the  top  of  the  master _denslieat.s  script: 


waittime 

nproc 

ntime 

nptile 

nqueue 

nproject 


wait  time  to  ping  queue  (default:  5  min) 

number  of  processors  for  Gaussian  calculation  (default:  32) 

queue  limit  for  Gaussian  calculation  (default:  4  h) 

number  of  cores  per  node  (default:  8  on  Harold) 

queue  to  use  (standard,  debug,  challenge,  background) 

project  account  code  (default:  determined  automatically) 


Please  remember  to  comment  out  any  duplicates  by  inserting  “##”  at  the  beginning  of  the  line. 

If  the  Gaussian  calculations  terminate  after  successful  submission  with  an  error  “died  without 
ever  signing  in,”  this  is  a  sign  that  the  secure  shell  is  corrupted  in  some  manner.  Perform  the 
following  in  order  to  obtain  a  successful  Gaussian  completion: 

Delete  the  $HOME/.ssh/  folder  with  the  command 

»  rm  -f  $HOME/.ssh/ 

and  wait  -10-30  min  for  it  to  re-appear.  Copy  the  following  commands  into  the  config  file  into 
the  new  $HOME/.ssh/  directory. 

Host  * 

ForwardXll  yes 
LogLevel=quiet 
Host  *  Banner  no 
StrictHostKeyChecking  no 
Host  * 
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ForwardXllTrusted  yes 
ForwardAgent  yes 

This  should  allow  successful  submission  of  future  script  calls  to  Gaussian. 


3.  Results  and  Discussion 


Upon  submission,  the  script  will  being  outputting  information  to  the  screen  (unless  the  standard 
output  is  piped  to  a  temporary  file  through  the  >  command  (as  denoted  previously)). 

Additionally,  up  to  three  output  files  will  be  generated  (along  with  a  host  of  intermediary  files). 
These  three  output  files  are  the  informational  output  file  ( arl_clensheat.err ),  the  heat  of  formation 
output  file  ( heat_of ^formation. out ),  and  the  density  output  file  (density. out). 

The  arl_densheat.err  file  is  a  general  message  file,  containing  what  has  been  done  along  the 
course  of  the  script  (e.g.,  files  built,  submitting  of  Gaussian  jobs  with  associated  queue 
information,  construction  of  ED  AT  files,  etc.).  If  the  script  fails  for  any  reason,  the 
arl_densheat.err  file  should  be  checked  first,  as  it  will  contain  information  such  as  whether  data 
was  not  copied  correctly  or  if  there  were  errors  in  the  Gaussian  calculations.  Dependent  on  the 
nature  of  the  error,  suggestions  may  be  offered  in  this  file. 

The  heatjof ^formation.  out  file  lists  the  name  as  given  by  the,  filename  (from  the  filename  .dat 
input  file),  the  chemical  formula  as  determined  from  the  input  data,  and  the  heat(s)  of  formation 
in  kcal/mol.  For  neutral  molecules,  the  gas,  liquid,  and  solid  phase  heats  of  formation  are 
computed,  for  both  the  atom  (a)  and  group  (g)  based  additivity  methodologies.  These  are 
denoted  as  Hfg  (a),  Hfg  (g),  Hti  (a),  Hfi  (g),  HfS  (a),  Hfs  (g).  For  ionic  materials,  if  the  compound  is  a  1:1 
salt,  two  methods  have  been  constructed  to  determine  the  solid  phase  heats  of  formation:  the 
Jenkins  and  Gutowski  methods  (2).  As  stated  above,  both  are  volume-based  methods  for 
determination  of  the  lattice  energies,  with  the  Gutowski  method  specifically  refit  for  1:1  salts.  If 
the  compound  is  any  other  salt,  the  more  general  Jenkins  method  will  be  the  sole  solid  phase  heat 
of  formation  listed.  A  sample  h ea t_o f  Jo n nali o n . o u l  is  presented  here  (if  the  text  wraps  to 
another  line,  please  expand  the  width  of  the  xterm  window): 

Name  Chemical  formula  Htg  (al  Hfg  (g)  Hfl  (a)  Hfl  (g)  Hfs  (a)  Hfs  (g)  (kcal/mol) 

nitromethane  (C  H3  N  02)  -18.555  -18.579  -26.246  -26.270  -32.740  -32.764 

15-DA4MeT_n  ([C2  H7  N6]l+)  ([N  03]1-)  Ionic  H((s):  33.744  Ionic  Hf(s)  (Gutowski):  12.016 

tag2bta  ([C  H9  N6]l+)2  ([C2  H  N9]2-)  ([H2  0])2  Ionic  Hf(s):  163.398  (kcal/mol) 

The  density,  out  is  similar  to  the  he  at _of formation,  out  file,  with  the  name  and  chemical  formula 
listed,  followed  by  the  crystal  density  in  grams/cubic  centimeter.  The  output  will  be  in  this 
format,  regardless  of  whether  the  compound  is  a  neutral  or  ionic  material.  A  sample  density. out 
is  presented  as  follows: 
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Name  Chemical  formula  Density  (g/cc) 

nitromethane  (C  H3  N  02)  1.424 

15-DA4MeT_n  ([C2  H7  N6]l+)  ([N  03]1-)  1.608 

tag2bta  ([C  H9  N6]l+)2  ([C2  H  N9]2-)  ([H2  0])2  1.498 


4.  Conclusions 


Scripts  have  been  successfully  developed  to  ease  the  complex  determination  of  the  crystalline 
densities  and  heats  of  formation  for  both  neutral  and  ionic  systems.  Simple,  minimal  information 
input  files  are  the  sole  requirements  needed  to  generate  these  important  performance  predictive 
properties.  The  scripts  will  grant  the  non-expert  user  the  capability  to  predict,  a  priori ,  the 
crystal  density  and  heat  of  formation  before  allocating  time  and  resources  towards  synthesis.  To 
quote  a  beta  tester:  “Prior  to  this  script,  it  was  very  labor-intensive  to  be  able  to  get  the  numbers 
from  the  computer.  There  was  a  series  of  three  calculations  per  molecule  which  had  to  be 
submitted  only  after  the  previous  mn  was  complete  (due  to  issues  with  trying  to  get  Gaussian  to 
combine  the  jobs  this  was  the  work  around),  then  very  long  command  strings  and  specific  file 
formats  that  had  to  be  remembered  and  entered  correctly  and  in  the  correct  sequence  in  order  to 
get  reasonable  numbers.  It  would  not  have  been  possible  for  someone  that  had  not  been 
specifically  trained  in  using  the  system  and  already  had  a  reasonable  knowledge  of  command 
based  computer  use  (in  other  words  old  school  without  a  GUI)  to  use  the  previous  method  of 
calculating  the  heat  and  density. 

With  the  new  script,  the  process  is  infinitely  simplified  with  only  one  .dat  file  to  build  with  the 
atomic  co-ordinates  and  the  correct  functional  group  allocations,  and  the  script  will  build  all  of 
the  required  files,  submit  the  jobs,  and  do  the  after-calculation  number-crunching  without 
requiring  any  further  input  from  the  user. 

Not  only  does  the  script  drastically  simplify  the  process,  but  it  also  cuts  the  required  time  by 
more  than  half  and  all  but  eliminates  transcription  errors  that  could  have,  and  frequently  did 
occur  when  . . .  first  learning  the  old  system,  in  building  the  correct  files  for  the  correct  job  (22).” 
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Appendix.  Supplemental  Figures  and  Instructions  for  Use  of  kinit,  PuTTY , 
and  gaussview 


Several  steps  are  required  prior  to  running  the  ARL  scripts,  such  as  gaining  access  to  the  HPC 
machines,  logging  into  the  HPC  machines,  and  building  the  Cartesian  coordinates.  This 
appendix  provides  a  quick  guide  for  these  required  steps,  but  is  not  intended  to  be  exhaustive. 
Please  contact  the  local  help  desk  or  CCAC  for  assistance  if  further  guidance  is  required. 

Prior  to  logging  into  the  HPC  center,  a  Kerberos  ticket  is  required,  as  this  ensures  secured 
communication  between  the  local  computer  and  the  HPC  supercomputer  center.  Running  the 
kinit  software  will  open  a  window  similar  to  figure  A-l  (note  that  for  a  completely  new  session 
all  boxes  will  be  blank).  Under  the  “name”  box,  please  enter  the  log-in  name  provided  from  the 
Consolidated  Customer  Assistance  Center  (CCAC)  and  enter  under  the  “realm”  box  the  text  seen 
below  (as  of  April  2012).  Click  on  “Login”  and  follow  the  rest  of  the  instructions.  If  successful, 
a  green  ticket  should  appear  (as  seen  in  figure  A-l).  If  not,  please  contact  CCAC. 


Figure  A-l.  Kerberos  kinit  window  to  obtain  ticket  required  for  access  to  HPC  computers. 

Once  a  Kerberos  ticket  has  been  obtained,  an  xterm  window  is  required  to  connect  to  the  HPC 
center.  Open  the  PuTTY  software  (figure  A- 2)  and  on  the  “Host  Name”  box,  enter  the  Harold 
address  (harold.arl.hpc.mil”,  and  under  “saved  sessions”  enter  the  name  you  wish  to  call  that 
session,  then  hit  save.  Henceforth,  one  can  open  a  connect  to  that  node  either  by  double  clicking 
on  the  saved  session  or  clicking  on  the  saved  session,  clicking  on  “load,”  and,  finally,  clicking  on 
“open”  in  the  bottom  right. 
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Figure  A-2.  PuTTY  window  to  open  connection  to  HPC  computers. 

In  order  to  interface  with  a  graphical  user  interface  (GUI),  the  local  machine  must  be  hosting  an 
xserver  emulator.  If  not,  the  HPC  machines  will  not  be  able  to  transfer  the  graphics  successfully 
to  the  local  machine.  Two  common  xserver  emulators  are  Exceed3D  and  Cygwin-X.  Please 
contact  the  local  help  desk  for  installation  of  either  of  the  xserver  emulators  to  the  local  machine 
and  ensure  that  the  xserver  portion  of  the  software  is  installed  (for  example,  for  Cygwin-X,  the 
XWin  Server  is  the  required  portion  of  the  software).  Once  installed,  double-click  on  the  xserver 
and  leave  it  running  in  the  background  prior  to  using  any  HPC-hosted  GUI.  As  stated,  failure  to 
have  an  xserver  running  in  the  background  on  the  local  machine  will  cause  the  HPC-hosted  GUI 
to  fail. 

Once  access  to  the  Gaussian  program  package  has  been  obtained  (means  to  test  access  in  main 
text),  access  to  the  gaussview  program  is  also  granted.  This  tool  is  a  graphical  interface  for 
building  and  visualizing  molecular  systems.  The  online  reference  manual  can  be  found  here: 
http://www.gaussian.eom/g  tech/gv5ref/gv5ref  toc.htm;  and  the  brochure  can  be  found  here: 
http://www.gaussian.eom/g  prod/gv5b.htm.  To  run  gaussview  at  the  prompt  (denoted  “»”  in 
the  text),  type  the  following  commands: 

»gview 
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If  the  xserver  emulator  is  running  properly,  figure  A- 3  should  appear  on  the  local  machine 
screen.  Please  refer  to  the  reference  manual  for  an  in-depth  tutorial  and  guidance  on  usage. 


Figure  A-3.  Gaussview  main  window. 


For  construction  of  molecular  systems  (needed  to  obtain  the  xyz  Cartesian  coordinates  required 
for  the  .dat  input  file  for  the  script),  three  buttons  will  be  used  most  frequently:  the  element 
builder  (figure  A-4),  the  ring  builder  (figure  A-5),  and  the  R-group  builder  (figure  A-6). 
Clicking  on  a  button  will  open  the  appropriate  builder. 


Figure  A-4.  Button  to  select  for  access  to  the  element  fragment  builder  in  gaussview. 
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Figure  A-5.  Button  to  select  for  access  to  the  ring  fragment  builder  in  gaussview. 


Figure  A-6.  Button  to  select  for  access  to  the  R-Group  fragment  builder  in  gaussview. 

The  elemental  fragment  builder  (figure  A-7)  can  be  used  to  add  individual  atoms  or  change 
existent  atoms  into  different  ones.  Select  the  desired  atom  and  valency  (bonding  type),  and  then 
click  on  the  molecule  window.  If  the  mouse  is  placed  and  clicked  in  a  void  region,  the  selected 
atom  will  be  added  and  capped  with  hydrogens.  If  the  mouse  is  over  an  existent  atom  and  then 
clicked,  that  atom  will  be  replaced  by  the  current  selection  from  the  element  fragment  builder. 

Hint:  The  bonds  displayed  in  the  molecular  viewing  window  do  not  indicate  actual  bonds,  as 
they  are  rendered  solely  on  atom-atom  bond  distances.  Gaussview  has  no  knowledge  of  the 
electronic  structure  (and,  thus,  the  bonds)  and  is  only  displaying  an  estimate  of  the  bonding 
structure. 

Hint:  If  one  wants  to  merely  replace  an  existent  atom  with  another  without  changing  the  current 
bonding,  when  selecting  the  bonding  type  (bottom  of  the  element  fragment  builder),  select  the 
“Atom”  (far  left)  designation;  this  will  not  alter  any  existing  bonds. 
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Figure  A-7.  Element  fragment  builder  in  gaussview. 

Gaussview  has  included  a  library  of  common  ring  (figure  A-8)  and  functional  groups  (figure  A- 
9),  which  can  be  added  in  a  similar  manner  as  an  element.  Selecting  a  fragment  and  clicking  on 
a  void  space  in  the  molecule  window  will  add  that  fragment  and  terminate  any  required  bonds 
with  hydrogens.  Clicking  on  an  existing  atom  will  replace  that  atom  with  the  selected  group. 
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Figure  A-8.  Ring  fragment  builder  in  gaussview. 
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Figure  A-9.  R-Group  fragment  builder  in  gaussview. 
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Once  the  general  molecular  geometry  has  been  constructed,  modification  of  the  bond  lengths, 
angles,  and  dihedrals  is  possible  (figure  A- 10).  Once  selected,  a  pop-up  box  will  appear, 
allowing  the  relevant  atoms  or  groups  to  move  or  be  frozen. 

Hint:  A  common  issue  with  geometry  optimizations  (as  performed  by  the  script)  is  the 
occurrence  of  imaginary  frequencies  (i.e.,  maxima  on  the  potential  energy  surface  instead  of 
minima).  A  vast  majority  of  the  time,  this  can  be  avoided  by  assuring  that  there  are  no  angles 
equal  to  exactly  90°  and  no  dihedrals  equal  to  exactly  180°.  For  technical  reasons,  structures 
containing  these  angles  and  dihedrals  usually  do  not  perform  the  proper  gradient-directed  steps 
required  to  optimize  the  geometry  to  a  minima  structure.  Breaking  these  symmetries  very 
slightly  (e.g.,  90.05°  versus  90°)  is  all  that  is  required  in  order  to  obtain  a  reasonable  optimized 
structure  that  still  has  the  desired  structure  (to  within  the  accuracy  of  the  theory). 


Modify 

angle 


Delete 

atom 


Modify 

bond 


Modify 

dihedral 


Figure  A- 10.  Options  to  modify  structures. 

Once  the  final  desired  structure  has  been  constructed,  it  is  necessary  to  save  the  geometry  for 
transcription  to  the  .dat  script  input  file.  Figure  A-l  1  is  the  menu  displayed  when  saving  a  file  in 
gaussview.  It  is  advised  to  have  “Save  as:”  set  to  “Gaussian  Input  File”  versus  “Auto.”  Once  the 
“Save  as:”  has  been  changed,  be  certain  to  select  “Write  Cartesians”  in  the  lower  left.  Type  an 
appropriate  file  name  and  hit  “Save.” 
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A-  Save  Input  File 

Look  m:  _j/usr/people/ebyr d/temp/ 


□  formic  acid.com 


& 

0-0- 

0-0- 

0-0- 

iff 

jformic_acid.com 

Save 

Gaussian  Input  Files  ( * . g j f  *.com) 

Cancel 

Gaussian  Input  File 

Help 

Write  Cartesians  m  Append  Extra  Input  J  New  Molecule  Group 


/a 


Figure  A-l  1.  Gaussview  save  menu. 


Using  one  of  the  UNIX-supplied  text  editors  (e.g.  vi,  emacs ),  open  the  saved  .com  file  (example 
illustrated  in  figure  A- 12)  and  select  the  atomic  coordinates  and  paste  them  into  the  appropriate 
.dat  file. 
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§chk=adsa  chk 

t  hf/  -  g  geom=connectivity 
Title  Card  Required 


Figure  A-12.  Saved  Gaussian09  input  file. 

Hint:  Recall  that  for  ionic  compounds,  each  individual  subunit  of  the  complex  is  separate  from 
the  other.  Either  build  and  save  each  subunit  separately  (and  paste  into  the  .dat  file),  or  build  the 
entire  system  and  edit  the  .com  file  appropriately  (i.e.,  delete  all  data  in  .com  file  except  for  the 
atoms  and  coordinates  and  use  the  “View,  Labels”  menu  option  in  gaussview  to  group  atoms  into 
the  correct  molecular  subunit).  Please  refer  to  online  manuals  for  detailed  instructions  on  the  use 
of  vi  and/or  emacs. 


29 


Intentionally  Left  Blank. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


6-31G* 

6-31G** 


6-311++G(2df,2p) 


A 

v 

2 

a2 

P 

Pcrystal 

AH°f  (g) 
AH°f  (1) 
AH°f(s) 

AH  Lattice 

c 


C’ 


H 

N 


N’ 


Pople’s  Gaussian  double-zeta  polarized  basis  set  with  cl  polarization 
functions  on  each  of  the  atoms  Li  through  Ca 

Pople’s  Gaussian  double-zeta  polarized  basis  set  with  cl  polarization 
functions  on  each  of  the  atoms  Li  through  Ca  and  p  polarization  functions 
on  H  and  He 

Pople’s  Gaussian  triple-zeta  polarized  basis  set  with  2 cl  and /polarization 
functions  on  each  of  the  atoms  Li  through  Ca  and  2 p  polarization 
functions  on  H  and  He  along  with  s  and  p  diffuse  functions  on  each  of  the 
atoms  Li  through  F  and  s  diffuse  functions  on  H  and  He 

o 

Angstrom 
Balance  parameter 

Average  electrostatic  surface  potential 

Variance  of  electrostatic  surface  potential 

Electronic  density 

Crystalline  density 

gas  phase  heat  of  formation 

liquid  phase  heat  of  formation 

solid  phase  heat  of  formation 

Lattice  enthalpy 

Carbon 

Carbon  with  hybridization 

Hydrogen 

Nitrogen 

Nitrogen  with  hybridization 
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n3 

no2 

o 

O’ 

R2 

Upot 

V 

X 

ARL 

B3LYP 

CCAC 

CPU 

DFT 

EM 

ESP 

G3MP2(B  3LYP) 


g/cc 

GUI 

HPC 

kcal/mol 

KS-DFT 

rms 


Azide  group 
Nitro  group 
Oxygen 

Oxygen  with  hybridization 
Coefficient  of  determination 
Lattice  potential  energy 
Electrostatic  potential 
Unknown/unassigned  atom  type 
U.S.  Army  Research  Laboratory 

Becke  3  parameter  exchange  with  Lee- Yang-Parr  correlation  DFT 
functional 

Consolidated  Customer  Assistance  Center 
Central  processing  unit 
Density  functional  theory 
Energetic  Materials 
Electrostatic  surface  potential 

G3  variant  which  uses  MP2  instead  of  MP4  for  the  basis  set  extension 
corrections  and  uses  the  B3LYP  structures  and  frequencies  instead  of  the 
Hartree-Fock  values 

grams  per  cubic  centimeter 

graphical  user  interface 

high  performance  computing 

Kilocalories  per  mole  (unit  of  energy) 

Kohn-Sham  density  functional  theory 

root  mean  square  error 
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PICATINNY  ARSENAL  NJ  07806-5000 
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ATTN  TECHL  LIB 
WASHINGTON  DC  20375-5000 


35 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  NAVAL  AIR  WARFARE  CTR 

ATTN  CODE  474100D  A  MERRITT 
ATTN  CODE  474100D  M  MASON 
2400  E  PILOT  PLANT  RD  MS  5202 
CHINA  LAKE  CA  93555 

2  NAVAL  AIR  WARFARE  CTR 

ATTN  CODE  470000D  A  ATWOOD 
ATTN  CODE  470000D  S  BLASHILL 
2400  E  PILOT  PLANT  ROAD, 

STOP  5001 

CHINA  LAKE  CA  93555-6107 

2  NAVAL  AIR  WARFARE  CTR 

NAV  AIR  WEAPONS  DIV 
ATTN  CODE  4T4200D  R  CHAPMAN 
ATTN  CODE  4T4320D  T  PARR 
1  ADMINISTRATION  CIRCLE 
CHINA  LAKE  CA  93555-6100 

1  NAVAL  SURF  WARFARE  CTR 

ATTN  A  STERN 
3824  STRAUSS  AVE  STE  107 
INDIAN  HEAD  MD  20640 

1  NAVAL  SURFACE  WARFARE  CTR 
RSRCH,  DEV,  TEST,  AND  EVAL 
IHDIV 

ATTN  C  MICHIENZI 
4081  N  JACKSON  RD  RM  18 
INDIAN  HEAD  MD  20640-5116 

2  NSWC  DAHLGREN 

G33  PRECISION  AND  ADVANCED 
SYS  BR 

ATTN  L  STEELMAN 
ATTN  N  COOK 
6210  TISDALE  RD  STE  223 
DAHLGREN  VA  22448-5114 

1  NVL  SUR  WARFARE  CNTR 

ATTN  H  HAYDEN 
4081  N  JACKSON  RD  RM  18 
INDIAN  HEAD  MD  20640-5116 

1  OFC  OF  NAV  RSRCH 

ATTN  B  ALMQUIST 
875  N  RANDOLPH  ST  RM  653 
ARLINGTON  VA  22203-1927 


1  OFFICE  OF  NAVAL  RSCH 

ATTN  C  BEDFORD 
875  N  RANDOLPH  ST  RM  653 
ARLINGTON  VA  22203-1927 

1  AIR  FORCE  RSRCH  LAB 

ATTN  RWME  T  KRAWIETZ 
2306  PERIMETER  RD 
EGLIN  AFB  FL  32542 

1  AIR  FORCE  RSRCH  LAB 

ATTN  J  ZABINSKI 
AFRL/RXB  B654  RM  135  2941 
HOBSON  WAY 

WRIGHT  PATTERSON  AFB  OH 
45433-7750 

1  LAWRENCE  LIVERMORE  NATL  LAB 

PHYS  AND  LIFE  SCI  DIRCTRT 
ATTN  P  PAGORIA 
LIVERMORE  CA  94550 

1  LOS  ALAMOS  NATL  LAB 
ATTN  L  PERRY 

PO  BOX  1663  MAIL  STOP  C290 
LOS  ALAMOS  NM  87545 

2  LOS  ALAMOS  NATL  LAB 
DYN  ENERG  MATS  DIV 
ATTN  D  CHAVEZ 
ATTN  P  LEONARD 
LOS  ALAMOS  NM  87545 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLA  F  DE  LUCIA 
ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLB  J  GOTTFRIED 
BLDG  1101 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWM  P  PLOSTINS 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005 


36 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLA  C  MUNSON 
ABERDEEN  PROVING  GROUND  MD 
21005 

I  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLB  B  HOMAN 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005 

I I  US  ARMY  RSRCH  LAB 
ATTN  RDRLWMLB  B  RICE 
ATTN  RDRLWMLB  E  BYRD 
(10HCS) 

ABERDEEN  PROVING  GROUND  MD 
21005 

3  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLB  I  BATYREV 
ATTN  RDRLWMLB  J  MORRIS 
ATTN  RDRLWMLB 
R  PESCE-RODRIGUEZ 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLB  R  SAUSA 
ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLB  S  BUNTE 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 
ATTN  RDRLWMLB 
WD  MATTSON 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 
ATTN  RDRLWMLC 
G  SUTHERLAND 
BLDG  1185 

ABERDEEN  PROVING  GROUND  MD 
21005 


1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLC  M  BISS 
ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLC  S  AUBERT 
BLDG  1119B 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLC  V  BOYLE 
BLDG  1185 

ABERDEEN  PROVING  GROUND  MD 
21005 

2  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLD  J  RITTER 
ATTN  RDRLWMLD  RA  BEYER 
BLDG  390 

ABERDEEN  PROVING  GROUND  MD 
21005 

2  US  ARMY  RSRCH  LAB 

ATTN  RDRLWMLE  P  WEINACHT 
ATTN  RDRLWMLG  W  DRYSDALE 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 
ATTN  RDRLWMLH 
M  FERMEN-COKER 
BLDG  309 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 
ATTN  RDRLWML  J  NEWILL 
BLDG  390 

ABERDEEN  PROVING  GROUND  MD 
21005 

2  US  ARMY  RSRCH  LAB 

ATTN  RDRLWML  MJ  ZOLTOSKI 
ATTN  RDRLWMM  J  ZABINSKI 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005 


37 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WMP  A  B  RINGERS 
BLDG  1116A 

ABERDEEN  PROVING  GROUND  MD 
21005 

2  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WMP  G  R  EHLERS 
ATTN  RDRL  WMP  PJ  BAKER 
BLDG  309 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  D  M  NUSCA 
BLDG  390 

ABERDEEN  PROVING  GROUND  MD 
21005-5056 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  B 
J  CIF.ZA  K  -JF.NKTNS 
BLDG  390 

ABERDEEN  PROVING  GROUND  MD 
21005-5059 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WM  B  FORCH 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  B  N  TRIVEDI 
BLDG  207 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  C  B  ROOS 
BLDG  1119B 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  C  KL  MCNESBY 
BLDG  1185 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 


1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  C  K  SPANGLER 
BLDG  1119B 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  C  T  PIEHLER 
BLDG  309 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  C  W  SHERRILL 
BLDG  1119B 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

2  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  F  D  HEPNER 
ATTN  RDRL  WML  F  D  LYON 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  F  E  BUKOWSKI 
BLDG  1119B 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  H  T  G  BROWN 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WMP  G  S  KUKUCK 
BLDG  309 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  WML  A  W  OBERLE 
BLDG  390 

ABERDEEN  PROVING  GROUND  MD 
21005-5067 
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NO.  OF 

COPIES  ORGANIZATION 


1  US  ARMY  RSRCH  LAB 

ATTN  RDRLROP  R  ANTHENIEN 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  US  ARMY  RSRCH  OFC 

ATTN  RDRL  ROP  C  J  PARKER 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  US  ARMY  RSRCH  OFC 

ATTN  RDRL  ROE  V  R  HARMON 
POBOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

4  US  ARMY  RSRCH  LAB 

ATTN  IMALHRA 
MAIL  &  RECORDS  MGMT 
ATTN  RDRL  CIO  LL  TECHL  LIB 
ATTN  RDRL  CIO  LT  TECHL  PUB 
ATTN  RDRL  D  V  WEISS 
ADELPHI  MD  20783-1197 
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Intentionally  Left  Blank. 
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